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Abstract 

We calculate the energies of ground and three low lying excited states of confined helium 

; 

^ atom centered in an impenetrable spherical box. We perform the calculation by employing 

. £h variational method with two-parameter variational forms for the correlated two-particle 
"c/3 ! 

K^iwave function. With just two variational parameters we get quite accurate results for both 

; 

i ^j ground and excited state energies. 
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I. INTRODUCTION 



Recently considerable attention has been focused on the study of spatially confined 
atoms and molecules 

HQS. 

The confined atomic and molecular systems show 
substantially different ground state and response properties as compared to their free 
counterparts. The main reason for the spatially confined models of atoms and ions to 
attract tremendous amount of attention is their applicability to several problems of 



physics and chemistry. For example, atoms trapped in cavities, zeolite channel 



Q0 



or encapsulated in hollow cages of carbon based nano-materials, such as endohedral 
fullerenes ((J z| and in nanobubbles formed around foreign objects in the environment 
of liquid helium and under high pressure in the walls of nuclear reactors are all 
relevant to the confined atom model. The models of confined atomic and molecular 
system have also found application in the investigation of effect of high pressure on 
the physical properties of atoms, ions and molecules The study of confined 

atoms also provides insight into various properties of the quantum nanostructures like 
quantum dots or artificial atoms jll|, lr2| [. The detail discussion on these applications 
are available in several review articles 

The first model of confined (compressed) hydrogen atom in an impenetrable spher- 
ical cavity was proposed by Michels et al. [9j to simulate the effect of pressure on 
hydrogen atom and this model was employed to study the variation of static dipole 
polarizability with the effective pressure acting on the surface. In this model the 
boundary condition that the wave function vanishes at r = r c ( where r c is the ra- 
dius of impenetrable spherical box) is imposed on the solution of the Schrodinger 
equation. Later various physical properties of confined hydrogen atom, such as the 
modification of their atomic orbitals, energy levels, the filling of electronic shells, 
inear and nonlinear polarizabilities have been reported in the literature ( see Ref. 
jj and references theirin). Besides hydrogen atom, effect of confinement by an im- 
penetrable as well as non- impenetrable spherical box, on many-electron atoms have 



also been considered 



13 



ipenetrable spnerical DC 



1J, [la, [16|, ll7|, Il8f, 19, 20]. Helium atom being the simplest 
many-electron system, the confined version of this atom provides a lucid way to study 
the effect of confinement on the electron correlation which arises due to the coulomb 
interaction between the two electrons. Most of the studies on the confined helium 
atom are devoted to the calculation of the ground state energies and some averages 
and their evolution with the size of the spherical box. In majority of these studies 
Raleigh-Ritz variational method was employed with modified Hylleraas-type wave 



function (Hylleraas-type wave functions multiplied with appropriate cut-off factor) 
fulfilling the confinement boundary condition mentioned above 14 , ^| . Besides 
variational method, self consistent Hartree-Fock configuration interaction 
and a quantum Monte Carlo (QMC) methods have also been used to study the 
properties of helium atom and several isoelectronic ions confined in an impenetrable 
spherical box. 

In this paper we report calculations of energies of the ground state and some low 
lying excited states of helium atom confined at the center of an impenetrable spherical 
box and study the variation of the energies with the size of the sphere. We note here 
that in comparison to the calculation of ground state energy very few studies on the 
evolution of excited state energies with the size exist in the literature (20, 21]. There- 
fore, main emphasis of the present paper is on the effect of confinement on some of 
the low lying excited states of compressed helium atom. The low lying excited states 
considered in this paper are 3 S(ls2s), 1 P(ls2p) and 3 P(ls2p). Recently, Patil and 
Varshni have calculated the the energies of the above mentioned excited states 
by replacing the electron-electron interaction by an effective screening of the nuclear 
charge. The screening factor is then determined by using an interpolation between 
the expressions for large and small values of the confining radius. In the present pa- 
per calculations are performed by employing more accurate variational method with 
two-parameter correlated wave functions for both ground and the excited states. The 
correlated wave functions explicitly take the effect of electron-electron interaction into 
account and consequently expected to yield accurate results. The variational forms 
for the wave functions we use in this paper are gen eralization of the correlated wave 



22, 



23, 



24L |25| for free two-electron 



functions proposed by Le Sech and co-workers 
atomic and molecular systems. The generalized wave functions for calculations of 
confined helium atom are constructed by multiplying the Le Sech type wave func- 
tions with appropriate cut-off factors so that the confinement boundary condition is 
satisfied. In addition to the confinement boundary condition these wave functions 
also fulfill both electron-nucleus and electron-electron cusp conditions. The cusp con- 
ditions arise due to the Coulomb interaction between the charged particles and the 

n 

true wave functions of many electron systems must satisfy these conditions 26]. At 
this point it is important to note that all previous variational calculations involved 
more than two variational parameters. We demonstrate that the two-parameter cal- 
culations performed in this paper yield quite accurate results and match well with 



the results of some accurate calculations already exist in the literatures. 

The remaining paper is organized in the following manner. In section II we describe 
the theoretical methods employed in this paper. The section III is devoted to the 
discussion of the results. The paper is concluded in section IV. 



II. METHOD OF CALCULATION 



The non-relativistic Schrodinger equation for confined two-electron helium-like sys- 
tems with nuclear charge Z can be written as ( in atomic units) 



-^Vi - ^V 2 + v N (r u r 2 ) + v c (r 1 , r 2 ) + v conf {r 1 , r 2 ) 
where is the nuclear potential 



^(ri,r 2 ) = Eif)(r 1 ,r 2 ) (I) 



Z Z 

v N (r 1 ,r 2 ) = , (2) 

n r 2 

Vc represents coulomb repulsion between the electrons 

v c (n,r 2 ) = — , (3) 

n 2 

and the confining potential v con f due to an impenetrable spherical box of radius r c is 
given by 

I \ I r l' r 2 < r c ... 

v^f (r 1 ,r 2 ) = < (4) 
^ oo n , r 2 > r c 

To solve the above Schrodinger equation for the both ground and excited states we 
employ Raleigh-Ritz variational approach by finding the stationary solutions of the 
following energy functional 

/ V*(ri, r 2) -f V? - \V\ + v N + v c + v conf ip(r u r 2 )dr 1 dr 2 



E 



(5) 



/ ^*(ri,r 2 )^(ri,r 2 )dr 1 dr 2 
In order to perform variational calculation we need to make a judicious choice for 
the ansatz of two-particle wave function. To this end we generalize the variational 



form for wave function ip(ri, r 2 ) proposed by Le Sech and co-workers 
The generalized variational form of the wave function, complying with the boundary 
condition imposed by the confining potential, employed in this paper is given by. 



V(rx, r 2 ) = 0(r l5 r 2 )ft(n, r 2 ) fl - d Vi _ r l 



(6) 



This form of the wave function is inspired by the concept of semi separability, in- 



troduced bv Pluvinaee 



by tne concept ol 

. y. hhs 



271. Following Ref. 22, 23l 24 25l d>{Y^ . To) is chosen in 



such a way that it only includes the wave function of the free electrons in the field of 

nucleus and f2(ri,r 2 ) contains the dependence of inter-electronic distance r i2 needed 

to represent the correlation arising due to coulomb interaction between the electrons. 

Moreover, the presence of cut-off factor ^1 — ^\ ^1 — £§J ensures the satisfaction 

of boundary condition imposed by the confining potential on the wave function. In 

addition to the satisfaction of confining boundary condition the above form of the 

cut-off factor also guarantees the fulfillment of the electron-nucleus cusp condition by 

the correlated wave function. The satisfaction of cusp condition by the variational 

wave function leads to improvement in the rate of convergence of Rayleigh-Ritz varia- 

I I 

tional calculations and also yields more accurate energies [26]. We note here that the 
calculations performed with linear form of the cut-off factor ^1 — ^1 — ^ , lead- 
ing to trial wave functions that do not satisfy electron-nucleus cusp condition, give 
less accurate results for the ground as well as the excited state energies, specially in 
the strong confinement regime. Consequently, in this paper we employ the quadratic 
form of the cut-off factor as given by Eq. For detail discussion on the various 
asymptotic and cusp properties satisfied by the Le Sech form of the wave function 
we refer the reader to Ref. 25]. With this choice for the variational form of the cor- 
related wave function, the energy functional (0) reduces to a single multidimensional 
quadrature 



where fl = Q yl — £h — ^J-J and E denotes the energy of two-electron systems 
moving in the nuclear potential only. The forms of the functions <f>(ri, r 2 ) and Q(r\, r 2 ) 
for the ground and some excited states employed in this paper for calculating corre- 
sponding energies are presented in the next section along with the results they yield 
for confined helium atom. 

III. RESULTS AND DISCUSSION 

In this section we present the results of our calculations and compare them with 
the results already available in the literature. First we discuss the results for the 
ground state energy and its variation with the size of the spherical box followed by 
the results for some of the low lying excited states of confined helium atom. 






A. Ground state 



The form of the function 0(ri,r 2 ) for the ground state ^(ls 2 ) of helium atom is 
written as 

(f>{r 1} r 2 ) = A (ls(n)l S (r 2 )) Hl)/3(2) - a(2)/3(l)] (8) 

where A is the normalization constant, ls(r) represents the ground state hydrogenic 
orbital and a(i) and denote the spin- up and spin-down functions respectively. 
To perform the calculations we next choose the form of function tt(ri, r 2 , ?*i 2 ), which 
describes the correlation between the electrons arising due to the coulomb interaction 

as 

f2(ri, r 2 , ri 2 ) = cosh Ari cosh Ar 2 ^1 + ^r 12 e~ ari2 ^ (9) 

where A and a are the two variational parameters which are determined by mini- 
mization of the energy functional (JJJ). The above form of Q(r\, r 2 , ri 2 ) consists of two 
parts, namely, the screening part represented by the cosh hyperbolic functions and 
purely r 12 dependent correlated part. The r 12 part of the the function was proposed 
by Hischfelder |28| and its applicability has been demonstrated for several many- 



electron systems up to beryllium atom 



22 



23 



24, 



25, 29]. This function provides 



25] 



right description of electron-electron cusp condition. In accordance with Ref. 
we choose to represent the screening part of fi(r 1 ,r 2 ,r 12 ) by cosh hyperbolic func- 
tions which fulfill the cusp condition at the nucleus and also right behaviour at large 
electron-nucleus distances. However, in contrast to the Ref. we represent the 
screening part by product of cosh hyperbolic functions at r\ and r 2 instead of the 
sum of two functions. Although the product form of the screening part considered 
in this paper overestimates the screening, nonetheless we find that this form leads 
reasonably accurate result for the ground state energy of an uncompressed or a free 
helium atom. The difference between our result for the ground state of uncompressed 
helium atom and the one obtained with sum of two cosh hyperbolic functions [3| is 
of the order of 0.002 a.u.. Having assessed the accuracy of our ansatz for the ground 
state energy of free helium atom we now proceed with the calculation of energies for 
confined helium atom as a function of size of the confining spherical box. 

In Table I we present the results of our calculations for energies of 1 S'(ls 2 ) state 
of helium atom as a function of r c along with the corresponding results available in 
the literature. In order to check the accuracy of our results we make a comparison 
with the results of some very accurate calculations based on correlated wave function 



14, 



1^ and QMC Q| approaches. It can be clearly seen from Table I that our results 
are slightly lower than the corresponding numbers of Ref. 14% especially for small 
values of r c . It is important to note here that in Ref. (l4| ground state energies 
were obtained by employing three-parameter variational calculation with different 
form for the function of r 12 . On the other hand our results are slightly higher than 
the results obtained both with the thirteen-parameter (ten linear and three nonlinear) 
Hylleraas wave function variational calculation [19j and the most accurate QMC based 
method. In the strong confinement regime (that is for small values of r r ) the maximum 

nn 

difference between our results and those of Ref. [la, 1121] is of the order of 0.006 a.u.. 
The results presented in Table I clearly demonstrate that the two-parameter ansatz 
used in this paper to calculate the ground state energies of confined helium atom 
gives quite accurate results for wide range of size of the confining spherical box. 

B. Excited states 

Now we apply the variational approach discussed above to calculate the energies 
of some low lying singly excited states and study their variations with the size of 
the spherical box. Following, Ref. j3, Q] the 0(ri,r 2 ) part of the correlated wave 
function for the excited state Isnl (where n is the principal quantum and / is the 
orbital angular momentum quantum number of the state to which one electron is 
excited) is chosen as 

0(r 1; r 2 ) = A [ls(n)raZ(r 2 ) ± ls(r 2 )nZ(ri)] x*(h 2) (10) 

where nl(r) is hydrogenic orbital with quantum number n and I and x a (l,2) repre- 
sents the spinor part of the wave function. The spinor part can be easily constructed 
by using spin-up a(i) and spin-down (3{j) functions such that the total wave func- 
tion should be antisymmetric with respect to the interchange of spatial and spin 
co-ordinates of two electrons. For calculation of energies of excited states the screen- 
ing cum r 12 dependent part of the wave function is chosen as 

0(ri,r 2 ,ri 2 ) = (cosh Xr\ + cosh Ar 2 ) ( 1 + -ri 2 e~ ari2 J (11) 



2 

The screening part in the above equation (Eq. ((HI ) is motivated by the work of 
Ref. j^J ■ Note that the form of screening part of the wave function for excited state 
calculation is different from that of ground state calculation. We find that product 
form leads to less accurate results for the excited state of uncompressed or free he- 
lium atom. Conseauentlv. for excited state calculation we emolov the above form 



(Eq. (Ilip) which has already been shown to yield accurate results for the excited 

r~i 

states of two-electron systems [23j. By using 4>(ti, r 2 ) given by Eq. (fTUj) we calculate 
energies of some low lying excited states for which the values of quantum number 
I or different spin state automatically ensures the orthogonality These states are 
3 5'(ls2s), 1 P(ls2p) and 3 P(ls2p). Before presenting the results for the confined he- 
lium atom we note that with Eqs. (fTUj) and (JTTJ) we get quite accurate results for the 
excited states of uncompressed helium atom. For excited state 3 iS(ls2s) of uncom- 
pressed helium atom we get E = —2.1743 a.u. in comparison to the most accurate 
value of E — —2.1752 a.u. [3] (here we quote up to fourth-decimal place only). On 
the other hand, our calculation for l P(ls2p) and 3 P(ls2p) states of uncompressed 
helium atom give E = —2.1227 a.u. and E = —2.1281 a.u. respectively, whereas the 
corresponding values of energies from accurate calculations are E = —2.1238 a.u. and 
E = —2.1332 a.u. jsi | (here we quote up to fourth-decimal place only) respectively. 
The results for energies of three excited states mentioned above and their variations 
with respect to r c are presented in Table II along with the corresponding results from 
Ref. |20j| in parenthesis. The comparison of two results clearly shows that our num- 
bers for 3 S(ls2s) states are lower than the corresponding number of Ref. jstj for 
almost whole range of r c considered in this paper except for r c = 1.0 a.u.. Moreover, 
for 3 S'(ls2s) state two results differ significantly in the range r c = 2.0 — 9.0a. -a.. In 
contrast to 3 S(ls2s) state case, our results for 1 P(ls2p) state are slightly lower than 
the corresponding numbers of Ref. for all values of r c excepting r c = 1.0 a.u.. 
Finally, we note from Table II that unlike excited states 3 5'(ls2s) and 1 P(ls2p) our 
results for the excited state 3 P(ls2p) are lower than the corresponding results of Ref. 
I20I ] only in the range r c = 2.0 — 9.0 a.u.. Beyond r c = 9.0 a.u. results obtained by 
us are little higher than those of Ref. j^j]. The lower values of energies for excited 
states obtained by our variational calculation particularly in the strong confinement 
regime suggests that our results may be more accurate than those of Ref. j^J. 



IV. CONCLUSION 



In this paper we have calculated the energies for ground and some excited states of 
confined helium atom by employing Raleigh- Ritz variational method. To perform the 
variational calculation we have used two-parameter variational forms for two-electron 
correlated wave function taking into account the boundary condition imposed by 
the confinement. The results obtained by us for the ground state is quite accurate 



within the strong confinement region and match well with the corresponding results 
of other accurate calculations. For excites states not many results are available in the 
literature. We have made comparison of our results for the excited states with the 
numbers from an interpolation based approach published recently. The comparison of 
results shows that our variation method with two-parameter wave function is capable 
of giving quite accurate numbers for the excited states of confined helium atom for a 
wide range of values of confinement radius. We feel that our results for the excited 
states will be an important contribution to the spectroscopic properties of confined 
helium atom and accuracy of these results can be further tested by more sophisti- 
cated calculations. The study of confined two-electron systems can be extended in 
several directions. The calculations of response properties like linear and nonlinear 
polarizabilities will provide information about the interaction of confined system with 
external electromagnetic field. Currently we are studying the effect of confinement 
on the linear and nonlinear polarizabilities of confined helium atom. 
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TABLE I: Energies for the ground state 1 S(ls 2 ) of confined helium atom as a function 
spherical box radius. All numbers are in atomic units. 



r c 


Present 


Ref. [14] 


Ref. [19] 


Ref. [18] 


0.6 


13.3343 


- 


13.3183 


- 


0.7 


7.9320 


7.9491 


7.9255 


- 


0.8 


4.6157 


4.6225 


4.6106 


- 


0.9 


2.4670 


2.4691 


2.4633 


- 


1.0 


1.0183 


1.0186 


1.0159 


1.0142 


1.1 


0.0106 


0.0106 


0.0091 


- 


1.2 


-0.7079 


-0.7075 


-0.7087 


- 


1.3 


-1.2304 


-1.2295 


-1.2309 


- 


1.4 


-1.6167 


-1.6151 


-1.6172 


- 


1.5 


-1.9061 


-1.9040 


-1.9067 


1.9081 


1.6 


-2.1253 


-2.1229 


-2.1263 


- 


1.7 


-2.2928 


-2.2903 


-2.2944 


- 


1.8 


-2.4219 


-2.4193 


-2.4242 


- 


1.9 


-2.5219 


-2.5195 


-2.5249 


- 


2.0 


-2.5998 


-2.5977 


-2.6036 


-2.6051 


2.2 


-2.7088 


-2.7074 


-2.7141 


- 


2.4 


-2.7765 


-2.7760 


-2.7831 


- 


2.6 


-2.8191 


-2.8194 


-2.8266 


- 


2.8 


-2.8462 


-2.8472 


-2.8542 


- 


3.0 


-2.8636 


-2.8652 


-2.8718 


-2.8727 


3.5 


-2.8851 


- 


-2.8928 


-2.8935 


4.0 


-2.8931 


-2.8956 


-2.8997 


-2.9003 


4.5 


-2.8963 




-2.9020 




5.0 


2.8978 


-2.9004 


-2.9028 


- 2.9032 


5.5 


-2.8985 




-2.9031 




6.0 


-2.8990 




-2.9033 


-2.9035 


oo 


-2.8999 


-2.9024 


-2.9035 


-2.9037 



TABLE II: Energies for three low lying excited states of confined helium atom as a function 
of spherical box radius. Numbers in parenthesis are results of Ref. 2Clt j . All numbers are in 
atomic units. 





3 5(ls2s) 


1 P(ls2p) 


3 P(ls2p) 


1.0 


15.5451 


8.0312 


7.5265 




(15.050) 


(7.751) 


(7.680) 


2.0 


0.5862 


-0.3414 


-0.4907 




(0.9809) 


(-0.3334) 


(-0.3692) 


3.0 


-1.3679 


-1.5217 


-1.5758 




(-1.1193) 


(-1.5069) 


(-1.5312) 


4.0 


-1.8734 


-1.8598 


-1.8907 




(-17277) 


(-1.8499) 


(-1.8686) 


5.0 


-2.0473 


-1.9928 


-2.0126 




(-1.9615) 


(-1.9857) 


(-2.0012) 


6.0 


-2.1171 


-2.0547 


-2.0684 




(-2.0658) 


(-2.0489) 


(-2.0624) 


7.0 


-2.1477 


-2.0861 


-2.0964 




(-2.1166) 


(-2.0812) 


(-2.0934) 


8.0 


-2.1617 


-2.1029 


-2.1111 




(-2.1429) 


(-2.0987) 


(-2.1101) 


9.0 


-2.1683 


-2.1121 


-2.1191 




(-2.1570) 


(-2.1087) 


(-2.1195) 


10.0 


-2.1714 


-2.1172 


-2.1234 




(-2.1647) 


(-2.1146) 


(-2.1249) 


11.0 


-2.1729 


-2.1199 


-2.1257 




(-2.1691) 


(-2.1181) 


(-2.1281) 


12.0 


-2.1736 


-2.1215 


-2.1269 




(-2.1716) 


(-2.1202) 


(-2.1300) 


13.0 


-2.1740 


-2.1223 


-2.1276 




(-2.1731) 


(-2.1215) 


(-2.1312) 


14.0 


-2.1742 


-2.1227 


-2.1280 




(-2.1739) 


(-2.1223) 


(-2.1320) 



